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ABSTRACT 

We have modelled the rotational modulation of X-ray emission from T Tauri stars 
assuming that they have isothermal, magnetically confined coronae. By extrapolating 
surface magnetograms we find that T Tauri coronae are compact and clumpy, such 
that rotational modulation arises from X-ray emitting regions being eclipsed as the 
star rotates. Emitting regions are close to the stellar surface and inhomogeneously 
distributed about the star. However some regions of the stellar surface, which contain 
wind bearing open field lines, are dark in X-rays. From simulated X-ray light curves, 
obtained using stellar parameters from the Chandra Orion Ultradeep Project, we 
calculate X-ray periods and make comparisons with optically determined rotation 
periods. We find that X-ray periods are typically equal to, or are half of, the optical 
periods. Further, we find that X-ray periods are dependent upon the stellar inclination, 
but that the ratio of X-ray to optical period is independent of stellar mass and radius. 
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1 INTRODUCTION 

One of the recent results from the Chandra Orion Ultra- 
' deep Project (COUP) was the detection of significant ro- 
' tational modulation of X-ray emission from low mass pre- 
, main sequence stars. The detection of such modulation sug- 
gests that the coronae of T Tauri stars are compact and 
clumpy, with emitting regions that are inhomogeneously dis- 
tributed across the stellar surface, and confined within mag- 
netic structur es that do not extend out to much beyond a 
stellar radius (iFlaccomio et al.ll2005l) . There is also evidence 
for much larger magnetic loops, possibly due to the inter- 
action with a surrounding circumstellar disc (jFavata et al.l 
I2OO5I : ICiardino et alll2006r) . A model already exists for T 
Tauri coronae, where complex magnetic field structures con- 
tain X-ray emitting plasma close to the stellar surface, whilst 
larger magnetic loops and open field lines are able to ca rry 
accretion flows iJardine et al.ll2n06l : lOregorv et al.ll200(tl . 

In this paper we use surface magnetograms derived from 
Zeeman-Doppler imaging to extrapolate the coronae of T 
Tauri st ars using stellar param eters taken from the COUP 
dataset (jCetman et all 120051) . By considering isothermal 
coronae in hydrostatic equilibrium we obtain simulated X- 
ray light curves, for a range of stellar inclinations, and 
then X-ray periods using the Lomb Normalised Periodogram 
(LNP) method. We compare our results with those of 
IFlaccomio et al.l (120051) . who demonstrate that those COUP 
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stars which show clear evidence for rotationally modulated 
X-ray emission appear to have X-ray periods which are 
either equal to the optically determined rotation period 
{Px = Popt) or are half of it {Px = 0.5Pc 



opt) 



2 REALISTIC MAGNETIC FIELDS 

In order to model the coronae of T Tauri stars we need to as- 
sume something about the form of the magnetic field. Obser- 
vations suggest it is compact and inhomogeneous and may 
vary not only with time on each star, but also from one star 
to the next. To capture this behaviour, we use as examples 
the field structures of two different main sequence stars, LQ 
Hya and AB Dor determined using Zeeman-Doppler imag- 
ing iD onati fc^^ollie^Cameronl L99^ 1 19971 : 

lDonati,.199ft : „Donati et al.lll999l:lD^ati et al.ll20oir Bv us- 
ing both stars, we can assess the degree to which variations 
in the detailed structure of a star's corona may affect its 
X-ray emission. The m ethod for extrapolatin g the field fol- 
lows that e mployed by|j ardine et al.' ('20 02a|) and is further 
detailed bv lJardine et al.> (|200S|) and iGregorv et al.1 (|2006l) : 
thus we provide only an outline here. Assuming the mag- 
netic field B is potential then V x S = 0. This condition is 
satisfied by writing the field in terms of a scalar flux function 
^, such that B = — V^. Thus in order to ensure that the 
field is divergence- free (V.B = 0), ^ must satisfy Laplace's 
equation, = 0; the solution of which is a linear combi- 

nation of spherical harmonics, 
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(a) (b) 

Figure 1. Closed coronal structures showing field lines which contain X-ray emitting plasma, extrapolated from surface magnetograms 
of (a) LQ Hya and (b) AB Dor, for COUP source number 6 (M* = O.23M0,i?* = 1.6i?0,Prot = 9. Sid). The magnetic structures are 
compact and inhomogeneously distributed about the stellar surface. 



* = ^ ^ [ai^r + 6(„>r-('+^>] P,™(^)e'™^ (1) 

l — l m— — l 

where Pim denote the associated Legendre functions. The 
coefficients aim and him are determined from the radial 
field at the stellar surface obtained from Zeeman-Doppler 
maps and also by assuming that at some height Rs above 
the surface (known as the source surface) the field be- 
comes radial and hence Be{Rs) = 0, emulating the effect 
of the corona blowing; open fi eld lines to form a stellar wind 
jAltschuler <^ Newkirk| [l^^ . In order to extrapolate the 
field we used a modified v ersion o f a code o riginally devel- 
oped bv lvan BallegooiienT Cartled ge Priest (1998). 

For a given surface magnetogram we calculate the ex- 
tent of the closed corona for a particular set of stellar pa- 
rameters. This is achieved by calculating the hydrostatic 
pressure at each point along a field line loop. For an isother- 
mal corona and assuming that the plasma along the field is 
in hydrostatic equilibrium then, 

p = po exp j( ^sc^s^ , (2) 

where Cs is the isothermal sound speed and Qs the component 
of the effective gravity along the field line such that, Qs = 
g.B/\B\. The constant po is the gas pressure at a field line 
foot point and p the pressure at some point along the field 
line. The effective gravity in spherical coordinates for a star 
with rotation rate co is, 

g (r, 0, (j)) — ^— -\- (jj'^r sin^ 0^ cj^r sin cos 0, 0^ . (3) 

We can then calculate how the plasma /3, the ratio of gas 
to magnetic pressure, changes along each field line. If at 
any point along a field line f3 > 1 then we assume that the 
field line is blown open. This effect is incorporated into our 



model by setting the coronal (gas) pressure to zero when- 
ever it exceeds the magnetic pressure {f3 > 1). We also set 
the coronal pressure to zero for open field lines. The gas 
pressure, and therefore the plasma f3, is dependent upon 
the choice of po which is a free parameter of our model. 
If we assume that the base pressure is proportional to the 
magnetic pressure then po — KBq, a technique which has 
been used successfully to calculate mean coronal densities 
and X-ray emissio n measure s for the Sun and other main 
seque nce stars yardine et al. 2002a b) as well as T Tauri 
stars (jJardine et al.M2006r ). By varying the constant K we 
can raise or lower the overall gas pressure along field line 
loops. If the value of K is large, many field lines would be 
blown open and the corona would be compact, whilst if the 
value of K is small, then the magnetic field is able to contain 
more of the coronal gas. The extent of the corona therefore 
depends both on the value of K and also on Bq which is 
determined directly from surface magnetograms. For an ob- 
served surface magnetogram the base magnetic pressure Bq 
varies across the stellar surface, and as such so does the 
base pressure po at field line foot points. By con sidering 
stars from the COUP dataset Ijardine et al.l (|200(t) obtain 
the value of K which results in the best fit to observed X- 
ray emission measures, for a given surface magnetogram. We 
have adopted the same values in this paper. We note that we 
make a conservative estimate of the location of the source 
surface by calculating the largest radial distance at which 
a dipole field line would remain closed, with the same av- 
erage field strength. Fig. Q shows examples of the complex 
magnetic field geometries which we consider. 

The coronal field geometries which we consider in this 
work have been extrapolated from surface magnetograms of 
the young main sequence stars AB Dor and LQ Hya. In fu- 
ture it will be possible to use real T Tauri magnetograms 
dervied from Zeeman-Doppler images obtained using the 
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ESPaDOnS instrument at the Canada-France-Hawaii tele- 
scope. However, in the meantime, the example field geome- 
tries shown in Fig. Q] capt ure three essential features of 
T Tauri coronae. First, Jardine et al.l (|200(tl have already 
demonstrated that the field structures which we consider 
here yield X-ray emission measures and mean coronal den- 
sities which are consistent with values obtained during the 
COUP - the largest avai lable dataset of X -ray properties of 
young stellar objects ( G etman et al-lEoOST) . Second, the fact 
that rotational modulation of X-ray emission was at all de- 
tected during the COUP automatically led to the conclusion 
that the dominant X-ray emitting regions must be compact 
and unevenly distributed across the star (|Flaccomio et al.l 
120051 ). It can be seen from Fig. Qthat the X-ray emitting 
closed coronal field lines do not extent out to much beyond 
the stellar surface (the coronae are compact) and that the 
emitting regions are inhomogeneously distributed about the 
star (the coronae are clumpy) , in immediate agreement with 
the COUP results. Further, as we discuss below, the field 
geometries considered in this paper do give rise to signifi- 
cant rotational modulation of X-ray emission. It can also be 
seen from Fig. Q] that some regions of the stellar surface do 
not contain regions of closed field, and would therefore be 
dark in X-rays. Such regions contain wind bearing open field 
lin es. Third, observations of var ious classical T Tauri stars 
bv lValenti Johns- Krud (120041) show the line-of-sight (lon- 
gitudinal) magnetic field components measured using photo- 
spheric absorption lines, are often consistent with a net cir- 
cular polarisation signal of zero. However, strong fields are 
detected from Zeeman broadening measurements. Such po- 
larisation measurements trace the surface field structure of T 
Tauri stars and immediately imply that such stars have com- 
plex and highly structured coronae. If the surface of the star 
is covered in many regions of opposite polarity this would 
give rise to a net polarisation signature of zero, with con- 
tributions to the overall signature from regions of opposite 
polarity cancelling out. Therefore the surface field must be 
highly complex and multipolar in nature, with many closed 
loops confining X-ray emitting plasma close to the surface of 
the star. However, it is also worth asking whether the coronal 
structure of classical T Tauri stars are likely to be funda- 
mentally different from those of weak line T Tauri stars, a 
point which we discuss in §4. 



Although we cannot be certain whether or not the 
magnetic field structures extrapolated from surface magne- 
tograms of young main sequence stars in Fig. Q] do repre- 
sent the magnetically confined coronae of T Tauri stars, 
they do satisfy the currently available observational con- 
straints. They reproduce X-ray emission measures and coro- 
nal densities which are typical of T Tauri stars, and the 
field structures are complex, confining plasma within un- 
evenly distributed magnetic structures close to the stellar 
surface. However, as such structures only represent a snap- 
shot in time of the coronal field geometry, we have chosen 
to consider two different field topologies. This allows us to 
model how much of an effect the field geometry has on the 
amplitude of modulation of X-ray emission and on resulting 
X-ray periods. 



3 X-RAY PERIODS 
3.1 X-ray light curves 

Once we have determined the structure of the closed corona 
we calculate the X-ray emission measure (EM) as described 
in Appendix A, and plot the ratio of EM to the maximum 
EM against rotational phase. Fig. |21 shows such plots using 
both the LQ Hya and AB Dor surface maps, for three incli- 
nations. The corresponding X-ray images are shown in Figs. 
13 and 31 from which it can be seen that emitting regions 
are compact and unevenly distributed across the stellar sur- 
face. However, some regions (which would contain open field 
lin es) are dark in X - rays, in agreement with the conclusions 
of iFlaccomio et all (|2005f) . that the saturation of activity 
cannot be due to the filling of the stellar surface with active 
regions. 

For the LQ Hya-Uke coronal structure, although there 
are emitting regions across the stellar surface, there are two 
dominant emitting regions in opposite hemispheres, one of 
which is brighter in X-rays than the other. The effect of this 
is most apparent for large inclinations. For an inclination 
oi i — 30° we see one minimum in the X-ray "light curve" 
(see Fig. |2(a)J, however as i is increased to 60° we see the 
emergence of a second minimum, almost completely out of 
phase with the first. For i — 90° both of the dominant emit- 
ting regions go into eclipse as the star rotates, giving rise 
to the two distinct minima. As one of these dominant emit- 
ting regions is brighter in X-rays, however, the X-ray EM is 
further reduced when the brightest region goes into eclipse. 
The amplitude of the modulation is around 50%. For an in- 
clination of 0° there is no rotational modulation of X-ray 
emission, as expected. 

The AB Dor- like coronal structure is more complex and, 
in contrast to the LQ Hya-like coronal structure, there are 
many bright emitting regions (see Fig. m. Fig. |2(b)1 shows 
the resulting X-ray light curves, and in this case we see little 
rotational modulation of the X-ray emission for i — 90°, but 
see modulation of around 50% for smaller inclinations. It is 
already clear from the shapes of the X-ray light curves that 
the X-ray period depends on the stellar inclination. 



3.2 X-ray and optical periods 

W e consider the same 233 COUP sources used in the analysis 
of IFlaccomio et all (|2005l) . however, as our model requires 
estimates of the stellar parameters (mass, radius, rotation 
period and measurements from which we can infer a coronal 
temperature) we have omitted those sources without esti- 
mates of these parameters. This leaves a total of 183 COUP 
sources with a lower coronal temperature and 141 where a 
higher temperature component is also given in the COUP 
dataset. For each COUP source we have generated an X-ray 
light curve (see Fig.|H]for three examples) using both the LQ 
Hya and AB Dor surface magnetograms, and for inclinations 
of 30°, 60° and 90°. In the subsequent discussion we only 
consider the lower coronal temperature. Our results for the 
higher coronal temperature are similar and are discussed in 
§4. Each simulated light curve covers 13.2 days, to match the 
COUP observation time. COUP light curves also contain 5 
gaps where the Chandra satellite passed through the Van 
Allen belts. Therefore we have considered both continuous 
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Figure 2. The variation in X-ray emission measure (EM) with rotation phase for (a) the LQ Hya-hke and (b) the AB Dor-hke coronal 
structures as shown in Fig. Q for inchnations of 30° (dashed) , 60'^ (dotted) and 90'^ (solid). There is clear rotational modulation of X-ray 
emission. 



(a) (b) (c) (d) 

Figure 3. X-ray images obtained using the LQ Hya-like coronal structure from Fig. |l(a)] for an inclination 90° (shown as the solid line 
in Fig. |2(a)] l. Emitting regions are inhomogeneously distributed across the stellar surface and are confined within magnetic structures 
close to the star; some regions are, however, dark in X-rays. This particular surface map has two dominant emitting regions in opposite 
hemispheres, which at this inclination, go into eclipse as the star rotates giving rise to the two minima in Fig. |2(a)] (solid line). Image 
(a) is for a rotational phase of 0.1, where the brightest of the two dominant emitting regions is in eclipse, (b) 0.4, where both of the 
dominant regions are visible, (c) 0.65, where the brightest emitting region is visible and the other is in eclipse and (d) 0.85, where once 
again both of the dominant regions can be seen. 



(a) (b) (c) (d) 

Figure 4. As Fig. |3lbut using the AB Dor-like coronal structure from Fig. |l(b)] and for an inclination of 30° (shown as the dashed line 
in Fig. |2(b)] l. Image (a) is for a rotational phase of 0.15, (b) 0.35, (c) 0.6, where the amount of X-ray emission is at a maximum due to 
the brightest emitting region (towards the bottom of the image) being visible and (d) 0.85. The coronal structure is more clumpy than 
the LQ Hya-like field, with many bright emitting regions across the stellar surface. 
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Figure 5. Comparison between X-ray periods Px and optically determined periods Popt for stars in the COUP dataset using (a) the 
LQ Hya-like (b) the AB Dor-like coronal structures shown in Fig.Q for inclinations of 30° (black circles), 60° (blue triangles) and 90° (red 
crosses). The solid lines represent Px = [0.5, l,2]Popt- The ratio Px/Popt is independent of stellar mass and radius. 



light curves (see Fig. |6(a)| >, and light curves which have gaps 
at approximately the same observation times (and lasting 
for approximately the same durati on) as the gaps in COUP 
X-ray light curves (see Fig. |6(b)| >. We tried sampling our 
simulated X-ray light curves every 2000s, 5000s and 10,000s, 
to match the analysis of Flaccomio et al. (2005), but found 
that our results were similar in all three cases. The results 
presented here refer only to a time sampling of 10,000s. 

We have calculated X-ray periods using the Lomb Nor- 
malised Periodo gram (LNP) me t hod s o that our analysis 
matches that of iFlaccomio et all (120051) as closely as possi- 
ble. The observed X-ray light curves of COUP sources are 
contaminated with many flares, with a great variety of flar- 
ing behaviour, whereas our simulated light curves do not 
suffer from this problem. Therefore, when calculating X-ray 
periods, our false alarms probabilities (FAPs) indicating the 
significance of peaks found in the LNP power spectrum are 
vanishingly small, with -40 < logio[FAP(%)] < -10. 

Fig. |51is a plot of our calculated X-ray periods against 
optica lly determined rotation periods from Flaccomi o et alJ 
i2005l). which is co l lection of rotation periods taken mainly 
from iHerbst et all (|2002f) . We have limited our search to 
X-ray periods of Px < 13(i and have chosen the scale of 
the axes to a llow a direct comparison with similar plots in 
IFlaccomio et al. (2005). For the LQ Hya-like coronal struc- 
ture, at any given inclination, the ratio Px/Popt is indepen- 
dent of stellar mass and radius, except for stars with large 
Popt (see the discussion below). We find that Px = Popt 
for inclinat ions of 30° and 60° , but that Px — 0.5Popt for 
i — 90° . kFlaccomio et all (|2005l) provide a qualitative argu- 
ment that the X-ray period would equal the optical period 
if there was one dominant emitting region on the stellar 
surface. They also point out that if there were two domi- 
nant emitting regions in opposite hemispheres, separated by 



180° in longitude, then it should be expected that the X- 
ray period would be half of the optical period. Although we 
agree with this point of view, our LQ Hya-like field structure 
does have two dominant emitting regions in opposite hemi- 
spheres, which for high inclinations yields Px — 0.5Popt but 
at low inclinations gives Px — Popt- Therefore the quali- 
tative picture of having two dominant emitting regions in 
opposite hemispheres giving rise to Px — 0.5Popt should 
come with the (perhaps obvious) caveat that it depends on 
the inclination at which the star is viewed. 

Results for the AB Dor-like coronal structure are simi- 
lar (see Fig. |5(b)| . We again find that the ratio Px/Popt is 
independent of stellar mass and radius for the i = 30° and 
60° cases, provided that Popt < lOd. Stars with Popt > lOc?, 
and in particular those which have optical periods longer 
than the COUP observing time (13. 2d), often do not have 
X-ray periods of either 0.5 Popt or Popt- This is most notice- 
able for the AB Dor-like coronal structure in the 30° case 
(see Fig. |5(b)| , where stars with the longest Popt are found 
to have X-ray periods of 0.7 Popt- This is a result of only 
considering the X-ray light curves over a duration equal to 
the COUP observing time. By extending our simulated X- 
ray light curves to cover twice the duration of the COUP 
observing time (so that ~ 2 complete rotations of stars with 
the longest Popt are covered by the light curves), we found 
that all X-ray periods in the i = 30° case clustered around 
Px = Popt- For the i = 90° case we find that most stars 
have Px = 0.5Popt, but five stars have Px = Popt- This is 
due to there being two significant peaks in the LNP power 
spectrum (see Fig. |6(c)| > and, depending on which peak is 
the larger, the X-ray period is either 0.5Popt or Popt- How- 
ever, when the AB Dor-like coronal structure is viewed at 
an inclination of 90°, the amount of rotational modulation 
of X-ray emission is at its smallest (only 13%) of all the 
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Figure 6. Light curves and Lomb Normalised Periodograms corresponding to inclinations of 30° (dashed) , 60° (dotted) and 90° (solid), for 
the coronal structure in Fig. |l(b)] The left hand column are the simulated light curves without gaps. However, when gaps are introduced 
into the light curves (right hand column) we find that the most significant peak in the LNP power spectrum generates the same X-ray 
period, albeit with lower significance. There are two significant peaks for the i = 90° case, one corresponding to Popt and the other 
O.SPopt- The solid vertical line is the frequency corresponding to the optical period. 



field structures and inclinations which we have considered; 
this can be seen from the light curve in Fig. |2(b)| Therefore 
the i = 90° case for the AB Dor-like field is an interesting 
result. It shows that it is possible to have a highly struc- 
tured corona, with emitting regions which are close to the 
stellar surface, with regions which are dark in X-rays, but 
when such a field structure is viewed at an unfavourable in- 
clination, little rotational modulation of X-ray emission is 
detected, and it is difficult to recover an X-ray period. How- 
ever, there is significant modulation detected for the i — 30° 



and 60° cases with Px — Popt- Therefore the X-ray periods 
we have determined are similar for both the AB Dor and 
LQ Hya-like coronal fields, and are dependent on the stel- 
lar inclination. The AB Dor field has many bright emitting 
regions making it difficult to qualitatively understand the 
exact shape of the resulting X-ray light curve. In the next 
section we further investigate the effects of inclination and 
the magnetic field structure on X-ray periods. 

We have also considered how having gaps in our sim- 
ulated X-ray light curves would affect the resulting X-ray 
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period (see Fig.E)). We found that the values of Px obtained 
from hght curves with gaps were identical to those derived 
from continuous light curves, with the only difference being 
that the main peaks in the LNP power spectra were of lower 
significance. 



3.3 Effect of inclination and coronal structure 

In order to further investigate the effects of inclination and 
the structure of stellar coronae, we randomly assign an incli- 
nation of between 5° and 90° to each of the COUP sources 
before calculating an X-ray period. Inclinations of smaller 
than ^ 5° produce little or no rotational modulation of X- 
ray emission and so are not considered. The results are plot- 
ted in Fig. 13 for both the LQ Hya and AB Dor- like coronal 
structures. We again find that X-ray periods are either equal 
to Popt or O.SPopt, with a few exceptions. 

Our simulations apparently highlight a well defined 
lower limit for X-ray periods of Px — 0.5 Popt- It should 
be noted, however, that we are limited by the number of 
surface maps which we have available; it could be the case 
that another surface magnetogram may yield X-rays peri - 
ods which are lower that 0.5Popt- IPlaccomio et al.1 (|2005l ) 
find stars with reliably detected rotational modulation of 
X-ray emission with Px < 0.5Popt but argue that such pe- 
riod detections may be spurious. In order to test if there is 
a well defined lower limit of Px = 0.5 Popt we generated a 
model coronal structure which had three emitting regions 
randomly positioned around the stellar surface. In the case 
where all three emitting regions were evenly spaced in lon- 
gitude, and a fixed (low) latitude (e.g. the three emitting 
regions placed at equal increments about the star's equator) 
we found X-ray periods of Px = 0.333Popt- Therefore it is 
possible to have a coronal field where the dominant emitting 
regions are distributed such that Px < 0.5Popt- However, we 
find that such field structures are rare in comparison to those 
whi ch yield Px = fO. 5 IjPo^ t . 

IPlaccomio et al.l (|2005fl also find a few stars with evi- 
dence for rotationally modulated X-ray emission with X-ray 
periods which are larger than the optical period, however 
their further simulations appear to indicate that these pe- 
riod detections (in particular those with Px > 2Popt) are 
Ukely to be spurious. We could not create any model which 
led to a Px value of much larger than Popt- 

The X-ray period depends on stellar inclination (see 
Pig. |8(a)] ). At low inclinations for the LQ Hya field, Px — 
Popt- As the inclination is increased a second peak in the 
LNP power spectrum becomes stronger and eventually, for 
i ~ 75°, both peaks are of almost equal significance. Por 
i > 75°, Px = 0.5Popt, with the secondary peak having 
become the most significant. This occurs when both of the 
dominant emitting regions enter eclipse at different times as 
the star rotates. At low inclinations the amplitude of mod- 
ulation is small (see Pig. |8(b)| > for the LQ Hya field, but 
is almost constant for i > 25°. Overall we find less mod- 
ulation for the AB Dor field. This is consistent with our 
finding from §3.2. When the AB Dor field is viewed at both 
low (i < 30°) and high (i > 75°) inclinations there is little 
rotational modulation of X-ray emission. Thus, simply hav- 
ing compact emitting regions inhomogeneously distributed 
about the stellar surface is not enough to generate large 
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Figure 7. Comparison between our calculated X-ray periods and 
observed optical periods. X-ray periods have been calculated for 
COUP stars from lPlaccomio et al.l (l200fJ) with randomly assigned 
inclinations. The lines represent Px = [0.5, 1, 2] Popt , with data 
for the LQ Hya-like (black circles) and AB Dor-like coronal struc- 
tures (red crosses). There is little difference between the different 
coronal field geometries. 

modulation of X-ray emission - the amount of modulation 
depends strongly on the stellar inclination. 



4 SUMMARY AND DISCUSSION 

By extrapolating surface magnetograms of young main se- 
quence stars we have shown that magnetic fields with a 
reaUstic degree of complexity reproduce the observed ro- 
tati onal modulation of X-ray emission from T Tauri stars 
(Placcomio et al.ll2005l) . This model has already been used 
to correctly predict X-ray EMs, mean coronal densities 
and the observed in crease in X-ray EM with stellar mass 
(jjardine et al.ll2QQ6^ . We find that X-ray emitting regions 
are inhomogeneously distributed around the stellar surface 
and ar e typically compa c t (< R ^). This agrees with the find- 
ings of IPlaccomio et all (|2005f) who point out that in order 
to explain the amplitude of the modulation of X-ray emis- 
sion which they detect, emitting regions must be compact 
in order to undergo eclipse. 

We found that the ratio Px / Popt is independent of stel- 
lar mass and radius, but that X-ray periods are dependent 
on stellar inclination. Por example, in the case where there 
are two dominant X-ray emitting regions in opposite hemi- 
spheres, we found that the X-ray period could be either 
0.5 Popt or Popt depending on how the star is inclined to 
the observer. In order words, the presence of two bright 
emitting regions in opposite hemispheres does not immedi- 
ately imply that Px — 0.5Popt- Purther, some coronal fields 
with compact emitting regions produce only small amplitude 
modulation of X-ray emission when viewed at unfavourable 
inclinations. 

We do not find X-ray periods of Px > "^Popt, unlike 
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Figure 8. (a) The variation in X-ray period as a fraction of the optical period with stellar inclination for the LQ Hya (black circles) 
and AB Dor (red crosses) coronal fields. Data represents the same stars as in Fig. Qwith the same randomly selected inclinations, with 
(b) the amplitude of modulation of X-ray emission. 



iFlaccomio et al.l (|20Q5fl . However, they argue that those pe- 
riod detections are likely to spurious, given the arduous task 
of removing flares from observed COUP X-ray light curves. 
Also the apparent lower limit of Px = 0.5Popt may not 
be of true physical origin. We have demonstrated this by 
constructing a corona which resulted in Px = 0.333Popt; 
however this sort of field structure is perhaps uncommon. 
Instead the majority of X-ray periods cluster around Px — 
[0.5, l]Popt- We found that Px (for a fixed inclination) was 
the same for both the LQ Hya and AB Dor-like coronal 
fields, despite the distribution of emitting regions being dif- 
ferent for the two cases. This suggests that for complex 
magnetic fields, where there are many emitting regions dis- 
tributed across the stellar surface, it becomes difficult to 
disentangle the contribution to the X-ray emission from any 
particular emitting region. Consequently, even though dif- 
ferent stars are likely to have different coronal structures, 
the X-ray periods of T Tauri stars appear to cluster around 
Popt and 0.5Popt- In other words, although different coronal 
structures give rise to different amounts of modulation at a 
fixed inclination, the X-ray period is likely to be the same. 

Throughout we have only considered the lower coronal 
temperatures derived from the COUP dataset. If instead we 
had used the higher coronal temperatures our results and 
conclusions remain the same. For a higher coronal temper- 
ature more field lines are unable to contain the coronal gas 
and are blown open to form a stellar wind, meaning that 
coronae are (slightly) more compact. This leads to small, 
but not important, differences in the resulting X-ray peri- 
ods, however the amplitude of modulation is larger (typi- 
cally 60% compared with 50% for the lower coronal temper- 
atures). The amplitudes of modulation wh ich we determine 
compa re well with the observed values of IFlaccomio et all 
HooB) of 20-70%. 



Although we have not considered it in this paper, the 
coronal structure of T Tauri stars may evolve on timescales 
that are comparable with the COUP observing time. For 
example a dominant X-ray emitting region may be lost due 
to a large coronal mass ejection which would open up a 
previously-closed region of field. This may then change the 
resulting X-ray period. However, it is still interesting that 
even with two different field structures such as those in Fig. 
Q] the X-ray periods are still typically equal to the optical 
period or are half of the optical period. 

We have not considered how the presence of a circum- 
stellar disc and active accretion influences coronal struc- 
tures, and the possible resulting effects on the X-ray pe- 
riod/amplitude of modulation. A comprehensive and de- 
tailed discussion of how accretion processes may influence X- 
ray em ission in T Tauri stars is provided bv lPreibisch et al] 
(|2005i) . The growing consensus, as was confirmed by COUP, 
is that X-ray emission from T Tauri stars does not depend 
upon the presence of a disc, but is influenc ed by active ac- 
cretion. For example, iFei gelson erall J2002) found X -ray ac- 
tivity lev els were the same f or sta rs with and without discs^. 
However, IFlaccomio et al.l (|2003r> found a difference in the 
X-ray activity levels of accreting and non- accreting stars. 
The process of accretion therefore influences the amount 
of X-ray emission detected, with stars that show evidence 
for accretion having X-ray luminosities which are a fac- 
tor of 2-3 smaller than those stars without discs and those 



^ It is important to remember that the presence of a disc does not 
immediately imply that accretion is taking place. Indeed many T 
Tauri stars show excess infrared emission, indicating the presence 
of a disc, but lack strong Ho; or Call emission, characteristic of 
active accretion. 
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stars which show evid ence for discs, but not for accretion 
dPreibisch et al.ll200.^ . 

Accretion may therefore affect the coronal structure of 
T Tauri stars, which in turn may have imphcations for ro- 
tational modulation of X-ray emiss ion, and X-ray periods . 
It is important to note however that lFlaccomio et al.l (120051 ) 
have searched for correlations between observed modulation 
amplitudes and the A(I— K) excess, a disc indicator, and 
EW(Call) , an accretion indicator , but do not find any. Fur- 
thermore, [Placcomi^^i] (|2005f) also search for dependen- 
cies on the X-ray period detection fraction with the indica- 
tors of discs and accretion, but again find none. This sug- 
gests that there is no difference in the spatial distribution 
of X-ray emitting regions in both classical and weak line T 
Tauri stars (with the former having both a disc and active 
accretion, and the latter neither). However. iFlaccomio et all 
(120051) did find a suggestion (but with a low significance of 
only about 1%) that stars with Px — 0.5Popt are preferen- 
tiall y active accretors. 

Ijardine et al.l (|2006f) have shown that stars (typically 
of lower mass) which have coronae that would naturally ex- 
tend to beyond the corotation radius would have their outer 
corona stripped by the presence of a disc. In order to check 
if this would affect the resulting X-ray period we determined 
the coronal structure of all the stars considered in §3 assum- 
ing that they were surrounded by a disc. Any field line which 
passed through the disc at, or within the corotation radius 
was assumed to have the ability to carry an accretion flow 
and was therefore considered to be "mass-loaded" and set 
to be dark in X-rays - a pro cess which is discusse d in detail 
bv lJardine et al.1 (l200(tl and lGregorv et al.1 ^00^. We then 
calculated X-ray emission measures with rotational phase 
and determined the amplitude of modulation of X-ray emis- 
sion and X-ray periods, using the same methods as discussed 
in §3. We found no difference in the values of Px and very 
little difference in the amplitudes of modulation. We there- 
fore conclude, exactly as the observations suggest, that the 
presence of discs does not influence rotational modulation 
of X-ray emission. However, active accretion might. 

We have yet to consider how accretion could influence 
rotational modulation of X-ray emission. As already dis- 
cussed above stars which are actively accreting show lower 
levels of X-ray activity than those which are non- accretors. 
Again, the pos sible explanations for this are discussed by 
IPreibisch et al.1 (Eo05i) . with the most likely being that mag- 
netic reconnection events in the magnetosphere cannot heat 
the dense material within accretion columns to a high 
enough temperature to emit in X-rays. Some models suggest 
that accretion will distort magnetic field li nes giving rise to 
instabilities and reconnection events (e.g. iRomanova et all 
I2OO4I ) . The energy liberated in such reconnection events may 
not be able to heat the high density material in accretion 
columns to a sufficiently high temperature in order to emit in 
X-rays. Thus accretion columns rotating across the line-of- 
sight to the star may be responsible for the observed reduc- 
tion in X-ray emission in accreting T Tauri stars compared 
with the non- accretors. The effect of this on X-ray peri- 
ods and modulation amplitudes will be considered in future 
work. A detailed accretion flow model will have to be con- 
structed in order to account for the attenuation of coronal 
X-rays by accretion columns rotating across the line- of- sight. 
Also, inner disc warps may have a similar effect at high stel- 



lar inclinations. This could change the shape of the X-ray 
variability curves shown in Figs. and |S1 with a reduction 
in X-ray emission at particular rotational phases, depending 
on the accretion flow geometry. 

Our model may allow for the simultaneous conflne- 
ment of X-ray emitting plasma in closed loops close to the 
stellar surface, but also for other field lines to stretch out 
into the inner disc and carry accretion fiows. iGregorv et al.l 
(|2002b) find that often open field lines can carry accretion 
fiows, which may allow accretion to infiuence the mag- 
netic field structure on the large scale, with field lines be- 
ing "mass-loaded" with disc material, but leaving compact 
coronae relatively unaffected by the process of accretion. 
This model would allow for the existence of large extended 
magnetic structures, inf erred to exist from the detection 
of large fiaring loops by iFavata et al.l (|2005l) , with the co- 
existence of compact coronae. However, much quantitative 
work remains to be done here. Another new accretion model 
may offer insights into the connection between X-ray emis- 
sion and accretion in T Tauri stars. An MHD model by 
I von Rekowski Piskunovl (|2006l) suggests that rather than 
material being loaded onto field lines at the inner edge of the 
disc, it fiows onto the star only at times when, and into re- 
gions where, the accretion fiow is able to overcome the stellar 
wind. In such a way this model predicts that accretion only 
occurs into regions of low field stren gth. It is also worth not- 
ing the observational results of »Flaccomio et all (|2005l) that 
indicate that rotational modulation of X-ray emission is not 
preferentially detected in either accreting or non-accreting 
T Tauri stars and that therefore the spatial structure of the 
X-ray emitting plasma is the same in both cases. Hence, 
although accretion does affect X-ray emission in accreting 
stars, and may distort the structure of stellar coronae, it 
may not affect the existence of rotational modulation. 
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where tm = [2$c(r^ - 1)]^''^[1 - r/rmV^^. Assuming that 
the pressure is uniform over the stellar surface we obtain the 
following expression for the emission measure as a function 
of r: 



EM{r) = 




(r2 - 1)1/2 



(A7) 

1/2 



exp 



2$^(^i-l^+2$c(r'-l) 



dr. 



where we remind the reader that all distances r have been 
scaled to a stellar radius. 

This paper has been typeset from a TgX/ I^T^ file prepared 
by the author. 



APPENDIX A: DETERMINING X-RAY 
EMISSION MEASURES 

In deriving expressions for the variation of the emission mea- 
sure with radial distance from the star, we begin with (|2|. 
If we scale all distances to a stellar radius, we can write 



p{r,0) =p*(l,6>)exp 



+ $c sin^ 6'(r^ - 1) 



(Al) 



where is the pressure distribution across the stellar 

surface and 



GM/R, 
kBT/m 



(A2) 

are the surface ratios of gravitational and centrifugal ener- 
gies to the thermal energy. The emission measure EM = 
J n^dV is then 



EM 



{kBTY 

/ pj(l,6>)exp [-2$c cos^ 6>(r^ - 1)] sin( 

J 6m. 



(A4) 



2 

r exp 



2$, 



+ 2$c(r - 1) 



dr. 



With the substitution /i = cos^ such that /^m — COS Un 
(1 — r/rm)"*^^^, the ^-integral can be written as 

/ pj(l,6')exp [-2$c/i'(r' - 1)] d/2 
Jo 

which, on using the substitution t — [2$c(^^ — l)]"*^^^/^, gives 
1 



(A5) 



[2$c(r 



11/2 



dt 



(A6) 



